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Development of a Nonlinear Probability of Collision Tool 
for the Earth Observing System 


David P. McKinley 1 

The Earth Observing System (EOS) spacecraft Terra, Aqua, and Aura fly in constellation 
with several other spacecraft in 705-kilometer mean altitude sun-synchronous orbits. All 
three spacecraft are operated by the Earth Science Mission Operations (ESMO) Project at 
Goddard Space Flight Center (GSFC). In 2004, the ESMO project began assessing the 
probability of collision of the EOS spacecraft with other space objects. In addition to 
conjunctions with high relative velocities, the collision assessment method for the EOS 
spacecraft must address conjunctions with low relative velocities during potential 
collisions between constellation members. Probability of Collision algorithms that are 
based on assumptions of high relative velocities and linear relative trajectories are not 
suitable for these situations; therefore an algorithm for handling the nonlinear relative 
trajectories was developed. This paper describes this algorithm and presents results from 
its validation for operational use. 

The probability of collision is typically calculated by integrating a Gaussian probability 
distribution over the volume swept out by a sphere representing the size of the space 
objects involved in the conjunction. This sphere is defined as the Hard Body Radius. 
With the assumption of linear relative trajectories, this volume is a cylinder, which 
translates into simple limits of integration for the probability calculation. For the case of 
nonlinear relative trajectories, the volume becomes a complex geometry. However, with 
an appropriate choice of coordinate systems, the new algorithm breaks down the complex 
geometry into a series of simple cylinders that have simple limits of integration. This 
nonlinear algorithm will be discussed in detail in the paper. 

The nonlinear Probability of Collision algorithm was first verified by showing that, when 
used in high relative velocity cases, it yields similar answers to existing high relative 
velocity/linear relative trajectory algorithms. The comparison with the existing high 
velocity/linear theory will also be used to determine at what relative velocity the analysis 
should use the new nonlinear theory in place of the existing linear theory. 

The nonlinear algorithm was also compared to a known exact solution for the probability 
of collision between two objects when the relative motion is strictly circular and the error 
covariance is spherically symmetric. Figure 1 shows preliminary results from this 
comparison by plotting the probabilities calculated from the new algorithm and those 
from the exact solution versus the Hard Body Radius to Covariance ratio. These results 
show about 5% error when the Hard Body Radius is equal to one half the spherical 
covariance magnitude. 

The algorithm was then combined with a high fidelity orbit state and error covariance 
propagator into a useful tool for analyzing low relative velocity/nonlinear relative 
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trajectories. The high fidelity propagator is capable of using atmospheric drag, central 
body gravitational, solar radiation, and third body forces to provide accurate prediction of 
the relative trajectories and covariance evolution. The covariance propagator also 
includes a process noise model to ensure realistic evolutions of the error covariance. This 
paper will describe the integration of the nonlinear probability algorithm and the 
propagators into a useful collision assessment tool. Finally, a hypothetical case study 
involving a low relative velocity conjunction between members of the Earth Observation 
System constellation will be presented. 
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The Earth Observing System began routine conjunction threat assessment of its member 
missions starting in 2004. It was quickly determined that a capability was needed to analyze 
conjunctions having nonlinear relative motion. An algorithm was developed that reduces 
the complex nonlinear relative motion into small linear segments that are easy to evaluate. 
The algorithm shows good agreement to the existing linear theory when used to analyze 
linear relative trajectories and for simple nonlinear cases for which solutions are known. 
The algorithm has been implemented in a software tool that accurately propagates the 
spacecraft states and covariances. The resulting 3D Collision Probability Tool is now being 
used by EOS analysts to study conjunction events. 


I. Introduction 

Starting in 2004, a collision risk identification and mitigation strategy was implemented for the Earth Observing 
System (EOS) series of spacecraft flown by Goddard Space Flight Center (GSFC). The 3 EOS spacecraft, Terra, 
Aqua, and Aura, fly with several other spacecraft in two constellations in 700 km mean equatorial altitude, sun- 
synchronous orbits. One of the key areas of interest was the capability to analyze and assess the collision risk during 
potential low relative velocity encounters between these spacecraft, since certain failure scenarios had been 
identified for the constellations that could lead to close approaches within a constellation [1], Since the spacecraft in 
each constellation are in essentially the same orbital plane, it is reasonable to expect that interactions between the 
constellation members can be described as low relative velocity. Although existing tools compute a high relative 
velocity Probability of Collision (P c ), the limits of applicability for that calculation as velocity transitions to a low 
relative velocity case are not well understood or documented. It was determined that a tool was needed that was 
capable of calculating the probability for low relative velocity encounters, not only so that low relative velocity 
encounters could be analyzed, but also so that an understanding of the transition from low relative velocity 
encounters to high relative velocity encounters could be analyzed. This paper documents the development of such a 
tool. 

A great deal of work has been done on calculating P c for high relative velocity encounters [2, 3, 4], The problem of 
low relative velocity has been investigated by Chan [5] and Patera [6], Chan’s method uses an analytical expression 
for the relative trajectory in order to derive the P c . Patera’s method uses a numerical scheme to evaluate the 
trajectory in increments to calculate the P c . 

This paper describes the algorithm and tools that have been developed for EOS to handle the calculation of P c 
between two spacecraft whose relative trajectory is nonlinear. The algorithm developed is similar to that of Patera 
[6], although the reduction of the volume integral into a contour integral is bypassed because the limits of integration 
that arise from the algorithm are easily expressed in rectangular coordinates. The derivation of the algorithm is 
presented, and then it is shown that the new algorithm is equivalent to existing algorithms for high relative velocity 
cases. The use of the algorithm to evaluate the breakdown of the high relative velocity calculations is briefly 
examined. Finally, a realistic . case study is examined through implementation of the algorithm within the 
commercial off-the-shelf astrodynamics software package FreeFlyer®. 
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II. Background 


The general structure of the probabilistic collision problem is shown in Figure 1. An object of concern, called the 
Primary, is predicted to come into close proximity with another object, called the Secondary. Without loss of 
generality, the Secondary’s position is taken as the origin and the Primary is visualized as moving along the relative 
trajectory. The combined geometric extent of the Primary and Secondary is described by a sphere of radius R^,, 
referred to as the combined Hard Body Radius (HBR). This sphere is attached to the Primary and sweeps out a 
volume, V, during the time of interest. A 3-dimensional relative probability density function (PDF), C(t), which is 
defined as the sum of Primary’s and Secondary’s position covariance matrices, is attached to Secondary. Generally, 
C(t) possesses a time-dependent shape and orientation due to the time variation in the covariance matrices of the 
separate objects. 


C(t) 



Figure 1: Typical Encounter Geometry. 

The probability of collision, P c , is defined as the integral of the PDF over the swept-out volume V [3]: 
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in some convenient set of coordinates. The vector l'’ 0 = {x 0 y 0 z Q } is the offset of 


the origin of the covariance from the center of the coordinate system. When the coordinate system and the 
covariance are centered on a common origin, the vector i' 0 is zero. 


In general the integral in Equation (1) is difficult to evaluate because of the time dependencies arising from the time- 
varying nature of C(t), and from the fact that different portions of C(t) fall within the Hard Body Radius at different 
times as the Primary evolves along the relative trajectory. 

Computation of Equation (1) can be simplified on those occasions when the relative motion is nearly rectilinear; this 
approximation is valid when the relative velocity between the Primary and Secondary objects is high. This relative 
motion, called the linear or 2D case, allows the well-known reduction of Equation (1) to a two-dimensional integral 
[2,3,4], Figure 2 demonstrates this linear configuration. A conjunction frame is constructed centered on the 
Secondary. The X -axis of the frame points to the Primary at the time of closest approach (TCA). The j> -axis is 
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along the direction of the relative velocity, and the z -axis completes the right handed coordinate system. This 
construction allows for the long axis of the cylinder V to be in the direction of the y -axis of the conjunction frame. 

Certain assumptions are made when solving the linear high relative velocity case. First, it is assumed that the 
encounter time is sufficiently small so that the individual covariance matrices and the resulting PDF are constant 
over the time of interest. Second, the high relative velocity means that the volume V will be several times longer 
than the 3-sigma value of the covariance in the y -axis. This equates to V essentially extending from negative 
infinity to positive infinity relative to the size of the combined covariance along the y -axis of the conjunction 
frame. Therefore, the integral along the y -axis is taken to be unity, and the problem expressed in Equation (1) may 
then be reduced to a convenient integral of a two-dimensional Gaussian PDF over the area, A (Figure 2), formed by 
the intersection of V with the conjunction plane: 




[f< 




dXclZ 
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where 



and C* is the 2-dimensional covariance matrix that results from projecting C(t) into the 


conjunction plane. 


However, the high relative velocity assumptions used in the construction of Equation (2) make it impractical to use 
in cases where the relative motion is nonlinear. The correct procedure would be to fall back on Equation (1) for a 
nonlinear trajectory. However, the nonlinear trajectory would create a complex geometry for the volume V, making 
integration difficult. In the next section an algorithm is presented that makes use of a series of simple geometries 
that approximate a nonlinear volume V. 
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III. Calculating the Probability of Collision along Nonlinear Relative Trajectory 

An algorithm similar to that given in Reference 6 has been implemented to investigate conjunctions with nonlinear 
relative trajectories. The algorithm differs from Reference 6 in that a new coordinate system in which to perform 
the numerical integration is developed. This limits of integration in the problem are easily defined in this new 
coordinate system, particularly when the problem of segment overlap is considered, as will be shown shortly. 

To accommodate an arbitrary geometry, the volume V will be dissected into smaller, simpler geometries. Consider 
a sphere that defines the combined geometries of the Primary and Secondary objects with radius Rhb as shown in 
Figure 3. During the encounter the sphere will similarly sweep out a long, curved, cylindrical volume. This volume 
can be thought of as made up of a series of smaller linear cylindrical volumes, V;. The geometric centers of the 
edges of the smaller cylinders are defined by the relative position of the Primary object with respect to the 

Secondary object at the current time and at a time step hi the future: Rrel, and Rrel m , respectively. By making the 
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time step between the two vectors small enough, the successive small linear cylinders can be made to approximate 
an arbitrary geometry defined by the nonlinear trajectory of the Primary relative to the Secondary. 



Figure 3: Dissection of complex geometry into small simple geometries. 

The cylinders, v ; , can be used as the limits of integration in Equation (1). The total probability of collision for the 
encounter will then be the summation over all vj. An additional advantage to this method is that the time steps may 
be sufficiently small such that the covariance is essentially constant over a single volume V;. However, the 
covariance can change from its value at subsequent time steps. 


The small simple geometries are easily integrated over using Equation (1) by constructing a convenient coordinate 
system at each time step, which is shown in Figure 4. The V i direction is defined as the difference between the two 

relative positions ^ REL and R I{EL ( . The N t direction is defined as the normal to the plane containing 

P P A A A A 

R rel and R rm i • The 5,. vector completes the right-handed system. The {V i N i B i } system is centered 
on li, m such that its origin does not coincide with the center of combined PDF. 


Vi 



Figure 4: Definition of integration frame. 


If the segments Vj are left as right cylinders, there is the possibility of overlap between segments as shown in Figure 
5. This overlap leads to integrating over the overlap volume twice, while not integrating over the excluded volume. 
When the HBR is large compared to the error covariance, this overlap issue can lead to incorrect probabilities [7]. 
To mitigate the overlap, each end of the cylinder Vj is cut off at an angle defined by the angle the adjacent segments 
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make to the current segment. These are angles 0i and 0 2 shown in Figure 6. There are four possibilities for the 
exact shape of each segment as shown in Figure 6. The exact geometry can be determined by comparing the cross 

A A /V A A 

products V xV l _ ] and V M xV t with the vector N i as described in the algorithm section. 

This method only accounts for overlap that occurs from nonlinearities in the in-plane relative motion. If there is a 

A A 

significant out-of-plane nonlinearity, i.e. N t is significantly different from N M , then additional measures would 

need to be taken to eliminate this overlap. Patera [7] suggests one methodology for dealing with significant three- 
dimensional overlap, and this topic remains one of the key areas for future work. 



Figure 5. Overlap of cylindrical segments. 



Cj out of the page 
C 2 out of the page 


C l into the page 
C 2 into the page 


Figure 6. Each segment v ; is truncated at the end by the angle between the current segment and the adjacent 
segments. Four shapes are possible depending on the relative angles of the adjacent segments. 
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IV. 3D Algorithm 


An algorithm was developed based on the concepts described in Section III. The following steps take the Earth 
Centered Inertial (ECI) states and covariances for the Primary and Secondary objects and compute the probability of 
collision. The Primary and Secondary states and covariances are propagated to a particular time i, the integration 
frame is constructed, and the P c associated with that segment of the trajectory is computed. This process continues 

and the overall P c is the summation of the individual P„ . 

c Cj 

1. The initial states and covariances at time i are labeled: J1 K $. p . , J2K lt s . , J2K C Pi , J2K C s . . The subscript 

P represents state and covariance for the Primary object and S represents the secondary. Propagate the 
J2000 inertial states and covariance to time i+1 and the inertial states also to time i+2 to obtain 

J2K K J2K j? J2K r* J2K J2K K J 2K K 
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2. Consfruct the transformation from J2000 to the integration frame at time i: 
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(4) 

(5) 

( 6 ) 
.(7) 
( 8 ) 
(9) 

( 10 ) 


3. When the individual position uncertainty covariance matrices J2K C p and 


J2K 


C o . are uncorrelated the 

A l 


uncertainty in the relative position between the primary and secondary object is described by the 
summation of the two matrices. Construct the combined covariance, transform it to the integration frame, 
and center it on the secondary object. 


re 'C,. = M™ K ( J2K C P + J2K C Si )M™ K 


( 11 ) 


4. Transform the miss distance vector from J2000 to the integration frame: 


INT 


EL, ^J2K J \ REL , 


wINT J2K 


l 


( 12 ) 


5. Compute the angles between adjacent segments 6 l and 0 2 : 

9 X = cos- 1 (V t 
0 2 =cos~ 1 (V i+r V i ) 


(13) 

(14) 


7 

American Institute of Aeronautics and Astronautics 



6. Compute the vectors C x and c 2 , which are used to determine the direction the trajectory is turning from 
segment to segment: 



7. By comparing the dot products c x • N t and c 2 • N { , the direction the trajectory is turning is known, and the 
limits of integration can be defined. The variables of integration in the { V. N t B t } system are v, n, and b. 
The limits for the v integral are given as a function of b and B l and 0 2 : 


If Cj • N ; > 0 and c 2 ■ N i < 0 


n 

b tan(~) < v < 


11K L,.r" K L, I + f-tan© 


A 

2 


(17) 


If Cj • < 0 and c 2 ■ N; < 0 


• b tan(-y) < v < \\ J2K $ RELm ~ J2K $re L , || + b tan(^) 


(18) 


If C, • Nj> 0 and C 2 • N { > 0 
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If c, • N t < 0 and c 2 - N t > 0 
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1. Compute the limits of integration in the and B t directions: 


( 20 ) 


■ 4 HBR 2 -b 2 <n< 4 hBR 2 - b 2 
- HBR <b< HBR 


( 21 ) 

( 22 ) 


9. Integrate the Gaussian distribution over the segment v, to determine the incremental probability associated 


with step i, where v,-, n t and b t are the components of the negative of the relative position vector: 
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v- V,. 
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n ~ n t 

b-b, 


( 24 ) 


10. Save R t _ j = R t . Return to Step 3 and repeat for the next time step until probability or trajectory conditions 

dictate that the loop ends. Discussion of when to terminate the loop is left to the following sections in 
which examples are given. 

1 1 . Accumulate the probability: 




(25) 


V. Comparison with Two-Dimensional Integration 

The algorithm presented in Equations (3) through (25), which will he referred to as the 3D method, makes no 
assumptions that would limit its use in linear as well as nonlinear relative trajectories. Therefore, using a linear 
relative trajectory, the probability can be computed using both the 3D method and the 2D method described in 
Equation (2). Simulated primary and secondary trajectories with 97-minute periods were generated with the 
Secondary object 500 meters directly below the Primary object in the radial direction at TCA. The angle between 
the two orbit planes was adjusted from 0 to 90 degrees to vary the relative velocity. A constant combined spherical 
position error covariance of 2 km 2 was used, and the probability of collision was calculated. 

Table 1 shows the results of the comparison of the 2D and 3D methods. As the angle between the two orbital planes 
decreases, the relative velocity at TCA decreases. For relative velocities as low as 13 m/s, the percent difference 
between the 3D and 2D probabilities is less than 1 percent, indicating good agreement between the 2D and 3D 
methods. As the relative velocities fall below about 2 m/s, the percent error grows quickly with the 2D method 
often giving probabilities an order of magnitude too conservative. 


Orbit Plane 
Angle 
(degrees) 

Relative 
Velocity (m/s) 

Pc 3D 

Pc 2D 

Percent 

Difference 

90.00 

10672.000 

9.393687E-05 

9.393670E-05 

1.809729E-04 

45.00 

5775.000 

9.330270E-05 

9.393690E-05 

6.751356E-01 

10.00 

1315.000 

9. 409870 E-05 

9.393670E-05 

1.724566E-01 

1.00 

131.700 

9.41 071 3E-05 

9.393670E-05 

1.814288E-01 

0.10 

13.175 

9.394736E-05 

9.393670E-05 

1.134807E-02 

0.01 

1.318 

1.644402E-04 

9.393670E-05 

7.505425E+01 

0.00 

0.270 

1. 489750 E-04 

9.393670E-05 

5.859084E+01 


Table 1: Comparison of 2D and 3D methods for relative trajectories with various relative velocities. 


The breakdown of the assumptions of the 2D method in this case happens below 13 m/s. Figure 7 shows the 
accumulated probability for the 5775 m/s and 0.27 m/s relative velocity cases. In the 5775 m/s case, the total 
Probability of Collision is “accumulated” in under two seconds. In the 0.27 m/s relative velocity case it takes nearly 
five hours for the probability to “accumulate”. The long dwell time of the 0.27 m/s case demonstrates the necessity 
of using a nonlinear algorithm to assess this probability. 

A general statement on the breakdown of the two-dimensional assumptions should not be inferred from this simple 
example. The failure of the two-dimensional assumptions is possibly dependent on the covariance size, shape (most 
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covariances are actually elliptical, not circular), and orientation. Further analysis is required to folly assess the 
transition between the two- and three-dimensional cases. Clearly, however, a breakdown of the two-dimensional 
case will occur, and a three dimensional tool is required to evaluate the collision probability in these cases. 


Probability of Collision 

5775 m/s Relative Velocity 
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“■ 0.00004 
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0.00001 
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2 

Elapsed Seconds 


3 




0.000150 

0.000125 

0.000100 

0.000075 

0.000050 

0.000025 

- 0.000000 


0.27 m/s Relative Velocity 



0:00 


3:00 


6:00 


Elapsed Hours 


Figure 7. Accumulation of the Probability of Collision for low and high relative velocities using 3D method. 


The 0.27 m/s relative trajectory depicted in Figure 7 is nonlinear, but is also essentially a fly-by trajectory. In fly-by 
trajectories one object approaches the vicinity of the other, reaches a minimum miss distance, and then departs the 
vicinity. In fly-by cases the conditions under which the integration loop can be ter m inated are clear. When the 
distance between the objects is large compared to the combined covariance the contribution of the trajectory 
segment Vj to the overall P c is negligible. Therefore, after the objects have left the vicinity of one another the 
probability no longer accumulates and the P c for this fly-by event is determined. Figure 7 then should not be 
interpreted as a time history of the P c ; rather it shows the summation of Equation 25 as the relative trajectory 
progresses. The P c must be tied to an event or time span, in this case a fly-by. The correct interpretation of Figure 7 
is that the steady state value of 1 .489E-04 shown is the P c for the given fly-by event. 
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VI. Nonlinear Trajectory 


A solution for the Probability of Collision has been given in Reference 6 for the case of a spherical hard body of 
radius r in a circular relative trajectory about the Secondary object with a spherically symmetric covariance with 
uncertainty a. The hard body radius centered on the Primary object sweeps out a torus of radius R about the 
Secondary as shown in Figure 8. The truth probability for this case is given by Patera [6]: 


2 2 

= k* exp| 


f 


y 


( R 2 +HBR 2 ) 

2c? 2 


I 


R^HBR 2 


\dx . 


(26) 


Chan [5] also derives a closed-form approximate solution for this case by assuming relative motion described by the 
Clohessey- Wiltshire equations. This approximation is given by: 


P c = Jhre - (R2+HBR)/2a2 R - WiR l. 

2o 


(27) 


Equation (26) was numerically integrated and compared to results from Equation (27) and the new 3D algorithm. 
The simulation parameters were set to R = 1 km and <r = 1 km. A one-degree step size was used around the torus. 
In this case the 3D algorithm was terminated after the primary object completed 360 degrees, or one relative orbit. 
The resulting P c is therefore the P c for the event defined as one orbit of the primary object about the secondary. The 
resulting probabilities for various values of hard body radius to uncertainty ratios (HBR/a) are shown in Table 2, 
and the percent differences are shown in Figure 9. The results show good agreement between the 3D method and 
the existing methods, with the 3D method differing by less than 1% from the exact solution of Equation (26). 


Primary Sphere 
with radius = 
HBR 



Secondary 

Object 


Figure 8: Circular Relative Motion of Primary about the Secondary. 
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HBRto 
Sigma Ratio 

3D 

Patera (Eq. 26) 

Chan (Eq. 27) 

0.01 

7.51446E-05 

7.58412E-05 

7.601 35E-05 

0.05 

1.87901 E-03 

1.89865E-03 

1.89805E-03 

0.1 

7.55945E-03 

7.57328E-03 

7.56382E-03 

0.2 

2.96574E-02 

2.99541 E-02 

2.98048E-02 

0.3 

6.58325E-02 

6.61437E-02 

6.54051 E-02 

0.4 

1.141 05E-01 

1.14537E-01 

1.12276E-01 

0.5 

1.72380E-01 

1.73008E-01 

1.67712E-01 

0.6 

2.38647E-01 

2.39024E-01 

2.28582E-01 

0.7 

3.09401 E-01 

3.09771 E-01 

2.91546E-01 

0.8 

3.82735E-01 

3.82306E-01 

3.53280E-01 

0.9 

4.53239E-01 

4.53694E-01 

4.10685E-01 

1 

5.20980E-01 

5.21 154E-01 

4.61069E-01 


Table 2. Comparison of 3D method with Patera’s exact solution and Chan’s approximation. 
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Figure 9. Percent difference between P c calculated using 3D method, Patera’s exact solution, and Chan’s 

approximation. 


VII. Case Study 

The 3D algorithm presented above was incorporated into the Conjunction Assessment and Mitigation (CAM) Tool 
Suite System, which was designed and built to support conjunction risk assessment and mitigation for EOS. The 
CAM Tool Suite was developed using the commercial off-the-shelf astrodynamics software FreeFlyer® (FF) and 
MATLAB®. FF provides the means to propagate the states and covariances using a wide set of force model options. 
Additionally, FF provides a built-in interface to MATLAB®, which is used to perform the numerical integration to 
determine the P c . 

A simulation of a hypothetical interaction between constellation members was developed. In this scenario, the 
Secondary object is slowly flying above the Primary object, and the relative motion is coplanar. At TCA the 
Secondary is approximately 50 m above the Primary and the relative velocity is approximately 1 m/s for the entire 
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encounter. Table 3 lists the orbital elements at TCA, and Table 4 lists the force model parameters used. The 
nonlinear relative trajectory in the Radial - In Track plane can be seen in Figure 10. The same initial covariance 
was used for each object in the Radial - In Track - Cross Track coordinate system: 


3 

0 

0 

0 

-0.010989 

0 

0 

34 

0 

-1.181196 

0 

0 

0 

0 

16 

0 

0 

0 

0 

-1.181196 

0 

0.04 

0 

0 

-0.010989 

0 

0 

0 

0.004 

0 

0 

0 

0 

0 

0 

0.02 


The time history of the Primary and Secondary object’s covariance can be seen in Figures 11 and 12. Figure 13 
shows that the P c computed for this encounter using the 3D method is 0.38877. Although this relative trajectory has 
a significant nonlinear component, this is again a fly-by trajectory and Figure 13 shows the summation of Equation 
25. The P c reported is the steady state value and represents the P c for this fly-by event. 

To substantiate these results, a Monte Carlo simulation was also performed. The same states, covariances, and force 
models were used as inputs. For each Monte Carlo trial, the Primary and Secondary states were randomly sampled 
according to their covariance matrices. The sampled states were then propagated forward in time. If the range 
between the states went below the HBR of 20 m, a conjunction was said to have occurred. The number of 
conjunctions divided by the number of trials was the resulting P c . Figure 14 shows the results of the Monte Carlo 
simulation. The P c computed by Monte Carlo was between 0.125 and 0.15, which is of a comparable magnitude to 
the results given by the 3D algorithm. This discrepancy is still under investigation. 


Parameter 

Primary 

Secondary 

a, km 

7077 

7078.8966 

e 

0 

0.0002665 

i, degrees 

98 

98 

Q, degrees 

221 

221 

w, degrees 

90 

90 

v, degrees 

0 

0 


Table 3. Primary and Secondary orbital elements at TCA 


Parameter 

Primary 

Secondary 

Mass, kg 

3239 

897 

Geopotential 

EGM96 4x4 

EGM96 4x4 

Cd 

2.2 

2.2 

Drag Area 

47.95 

8.3 

Atmosphere 

Jacchia Roberts 

Jacchia Roberts 

3 ra Body 

Sun, Moon 

Sun, Moon 


Table 4. Force model properties for case study. 
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Radial vs. In Track 


4.0 



- 10.0 - 7.5 - 5.0 - 2.5 0.0 2.5 5.0 7.5 10.0 

In Track 

Figure 10. Nonlinear relative trajectory of the Secondary object relative to the Primary object. 
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Figure 11. Primary object Covariance time history. 

Secondary RIC Error Sigma 
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Figure 12. Secondary object Covariance time history. 


14 

American Institute of Aeronautics and Astronautics 



Probability of Collision 
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Figure 13. Probability of Collision time history. 


Probability of Collision Vs. Iteration 

convergence check 




VIII. Summary and Future Work 

An algorithm based on [6] has been implemented to calculate the probability of collision between spacecraft with 
nonlinear relative trajectories. The complex geometries that can be created in these encounters are broken down into 
simple geometries that are easily used as integration volumes. The resulting algorithm compares extremely well 
with existing 2D theory when analyzing linear relative trajectories. In addition, the new algorithm has also been 
shown to give comparable answers to a published reference solution for the case of a highly nonlinear circular 
relative trajectory and compares within an order of magnitude with Monte Carlo results. 

The 3D tool is currently being used to analyze conjunctions with EOS and its fellow constellation members. While 
most conjunctions are high velocity, using the 3D tool allows the EOS Conjunction Assessment Analysts to begin 
gaining insight into the transition between linear and nonlinear relative trajectories. Additionally, should any low 
velocity encounters between constellation members or between EOS and other space objects occur, the tool is 
available to assess the collision risk. 
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